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Abstract 

A semiclassical approach for calculating shell effects, that has been 
used in atomic and plasma physics, is applied to describe electronic su- 
pershells in metal clusters. Using the spherical jellium model we give 
the analytical expression for the oscillating part of the binding energy of 
electrons as an explicit sum of contributions from supershells with quan- 
tum numbers 2nr + I, Srir + I, 4nr + I, ... This expression is written in 
terms of the classical characteristics of the motion of an electron with the 
Fermi energy in a self-consistent potential. The conditions under which 
a new supershell appears and the relative contribution of this shell are 
studied as a function of the cluster size and form of the potential. Specific 
calculations are performed for a square well. 

1. The electronic structure of atomic clusters has been intensively studied ex- 
perimentally and theoretically in the last two decades (see, for example, reviews 
in Refs. and [^). One general feature of the experimental mass-abundance 
spectra of TV-atom clasters is the existence of magic numbers: Clusters with 
these numbers of atoms are produced more abundantly. As N enlarges, the 
amplitude of these variations lessens, then in still larger clusters increases once 
again, and so on, i.e., oscillations with beats occur. 

The theoretical reproduction of the beating patterns one can see, for exam- 
ple, in the self-consistent calculations of sodium clusters in Ref. (s), where the 
spherical jellium model in local density approximation was used. There are two 

kinds of periodicity in the density of states and in the oscillating part SEsh of 
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the binding energy as a function of a cluster radius R ^ Ne {Ne is a number 
of valence electrons in cluster): the primary shell periodic structure and beat 
mode with a higher-order period. The similar results were performed in Ref. Q 
using the Woods-Saxon potential and in calculations for various metals in Ref. 
using nearly self-consistent potentials. 

For small clusters the interpretation of the numerical results is obvious: the 
extremum energy cusps occur at Imax-sheW closing and represent the magic 



1 



numbers (Imax is a maximum orbital number I in the cluster) . However this law 
breaks down for iVg > 100. 

The theory, explaining quantitatively the beating pattern of the level density 
for spherical cavity by a superposition of the contributions from closed classical 
trajectories of electrons, has been developed in the fundumental work of Balian 
and Bloch The detailed numerical calculations in Ref. |^ using more com- 
plex spherical cluster potentials reproduce the similar oscillations, one-electron 
levels s{nr,l) of the higher angular momentum states bunching in supershells: 
e{nr, I) ~ e{nr + 1,1 — K) with pscudoquantum numbers Krir + I. Here is a 
radial quantum number, K — 2,3,4, .... The integer number K under classical 
treatment (see, for example in Ref. ) is equal to a ratio of the frequencies in 
the radial and angular motion of the corresponding closed orbit. Ai K = 2 
pscudoquantum number is the same as a principal quantum number and char- 
acterizes the pendulating orbit going through the origin, if = 3 and K = A 
corresponds to the triangular and square orbits, respectively. 

One can suppose from here that a periodic-orbit-expansion (or supershell- 
expansion), obtained for a spherical cavity, is a particular case of the more gen- 
eral expansion for spherical cluster potentials. This conclusion is supported in 
the recent paper Q where such a generalization was carried out using Woods- 
Saxon potential by expanding on a parameter a/ R (a is a surface width) around 
the known results for a spherical potential well. It is of interest to investigate 
analytically the origin of supershells and the mechanisms leading to their ap- 
pearance for a potential of arbitrary form. 

In the present paper it is shown that this problem can be solved by a semi- 
classical method for distinguishing shell effects, previously applied successfully 
in atomic [|j , ||l^ , and plasma ||ll| , ||l^ physics on the basis of the Thomas-Fermi 
(TF) model. Although the TF model and its conventional variants ETF with 
quantum and exchange corrections (for application to clusters, see, for exam- 
ple, Refs. and [Q) give only the average dependences of all quantities on 
the number of particles, a refinement of this model makes it possible to account 
for the shell structure of the electronic spectrum. This refinement is based on 
the use of the Bohr-Zommerfeld quantization conditions and on the possibillity 
of performing the sum over quantum numbers analytically, provided that the 
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semiclassity parameter, which for clusters is proportional to Ne , is small. 

2. We furthermore use the expression derived in our future extended pub- 
lication for the correction to the binding energy of electrons (the atomic units 
are used): 



Here (5n(r) is the correction to the electronic density, which because of other 
effects goes beyond the initial model approach (for example, TF or ETF model), 
jjL is the chemical potential in the initial model, the correction i5n(r) is assumed 
to be small and is calculated using the initial self-consistent potential. 

We are interested in the contribution of the shell correction to the electronic 
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density of states (5ns/i(r,A') or to the number of states: 



where for a cluster with the closing 'T' -shells 

iV(Ai) = 2^(2? + l)0(^-£„^,O, 

71 J-, I 



(2) 



(3) 



and in the semiclassical approach the energy levels £^^,1 are determined from 
the quantization condition 



drpsi (r) = TT rij, + 
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Here S,i = S,x and p,i{r) = ^2{e - U{r)) - {I + l/2f /r^ = ^pl{r) - AVr^ = 
PeA are, respectively, the classical radial action and momentum of an electron 
with energy e and orbital angular momentum /. Simple calculations using the 
Poisson formula to replace the sums over quantum numbers rir and I by integrals 
make it possible to rewrite expression (3) as 



k,s — — oo 



dXX sm{27rkiy^x) cos(27rsA) 



(5) 



Here v^x = Ssx/tt, and determines the border of the phase area of the 
classically allowed motion of an electron with an energy e: z^^Ae — 0. In Eq.(5) 
the term with k — s — gives the TF result Ntf{iJ') and according to Eq.(2) 
the sum (5) without this term is equal to the desired quantity SNghifJ-)- 

Let's note now, that the supposed approach is agreed with a general, known 
from the nuclear physics, concept jl6| of separating the total energy of finite 
sistem as a function of the system size into a smooth part and a fluctuating 
correction. But here for calculating the last we base the simple expression (1), 
using initial smooth self-consistent potential or its Woods-Saxon and another 
fitting. 

3. The integral limits and points A of stationary phase make the main 
contribution in the integral over A in Eq.(5). These points are determined from 
the relation: 



dX 



< A < A^ 



(6) 



s 

A fc' 

A function i'^a(A) decreases monotonically and for all potentials U{r) that are 
finite at the origin the slope of the corresponding curve at A = is the same 



[0, 



dX 

The value of the derivative at A = A^ 

dX 



(6a) 



(66) 
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depends strongly on the form of the potential. For an oscillator = 1/2, for 
a square well v'^j^ = 0, and for the Woods-Saxon potential the value of i^J^ varies 
with increasing N, vanishing in the limit of a very large number of atoms. 

On this basis it follows that the relation (6) distinguishes in the sum over k 
the leading terms 



k={2+j)s, 



du. 



1 



dX 

j ~ 0? 1; •'•jmaxi 3 max 



2 + J 
1-2 



(7) 



The terms k = 2s {j — 0) must be studied separately, since in this case the point 
of stationary phase A = is also the lower limit of integration. As a result we 

obtain 

_ (^1)'' J cos(27r.s2t/^o) cos(27rsA^) 



Here 



2tts 



(8) 



Substituting expression (8) into Eq.(l) and integrating by parts to separate 

the terms which are of leading order in the semiclassicity parameter give a 
semiclassical formula for the shell correction to the binding energy of electrons 
in a cluster 
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(-1)^ j sin(27rs2i^^o) 



A^ sin(27rsA^) 



2 ^ {nsf \ (a2z.^o/a/x) (0.5 - i^',) {dXJdn) 



4Vi.j.(-ip- sin 



2tts[{2 + j)v^-^^+X, 



(9) 



as an sum of contributions from supcrshcUs with quantum numbers Uj = Krir + 
K = 2 + j, the quantization on the Fermi level being substantial. In Eqs.(8) 
and (9) the symbol j numerates a kind of the electronic orbit: j = corresponds 
to the linear pendulating orbit, the terms j > 1 {K > 3) are connected with 
the planar regular polygons, K being the number of their vertices. The integer 
value s is equal to a number of periods that the electronic trajectory {j, s) 
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includes, so the sum over s is the trajectory length-expansion for a j-orbit 
(compare with the orbits {X,i^) in Ref. |^). 

4. The proposed simple method makes it possible for any spherical poten- 
tial to determine the period and amplitude of oscillations associated with each 
supershell and to estimate their relative role in beating occurence. The results 
of such an analysis for a square well potential: 



U{r) 



-2e_F,r < R 
0,r>R 



R = r,N, 
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are displayed in Fig.l. Let's note that even though we start with the expression 
(3) for a cluster with the closing "Z"- 
shells our results for Ne > 100 check 
well with the results of the complete 
calculations in Ref. H. Fig.l shows, 
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that the first period ( A'e < 7) of the 
beats is determined by terms with j ~ 
0, 1, 2, the term with j — contribut- 
ing little. The account for term with 
j = 3 {K — 5) is needed to describe 

1 /3 

well the second period (A^e' < 13). Add- 
ing a term with j = 4 {K = 6) is sufh- 
cient to describe the behavior of SEgh 
in the entire range under study. This 
means that the actual value jmax is less 
that determined from Eq.(7) and cor- 
responds to the filled states, for which 
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5. For spherical potential we give 
the analytical expression for the oscil- 
lating part of the binding energy of 
electrons as an explicit sum of contri- 
butions from the leading periodic or- 
bits (supershells). We show that a 
number of the leading orbits in the ex- 
pansion depends strongly on the form 
of the potential and is determined by 
the value of derivative (6b)(see Eq.(7)). 
The more soft potentials will be dealt 
in detail in our extended publication 
where among other things small clus- 
ters and finite temperatures will be an- 
alyzed. 
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a) 



Shell energy correction 
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Figure 1: 

SEsh as a function of Ne'''- The 
unit of energy is the Fermi energy ep- 
b)-f) - Comparasion of a truncated 
supershell-expansions. 



5 



References 

[1] W. A. de Heer, Rev.Mod.Phys.65, 612 (1993) 

[2] M. Brack, Rev.Mod.Phys.65, 677 (1993) 

[3] O. Genzken and M. Brack, Phys.Rev.Lett 76, 3286 (1991) 

[4] H. Nishioka, K. Hansen and B. R. Mottelson, Phys.Rev. B42, 9377 (1990) 

[5] K. Clemenger, Phys.Rev. B44, 12991 (1991) 

[6] R. Balian and C. Bloch, Ann.Phys. 69, 76 (1971) 

[7] S. Bj0rnholm, in "Nuclear Physics Concepts in Atomic Cluster Physics" 
(Springer, 1992), pp.26-38 

[8] E. Koch, Phys.Rev. B58, 2329 (1998) ( [cond-mat /9803309| ) 

[9] D. A. Kirzhnits, Yu. E.Lozovik and G. V. Shpatakovskaya, Usp. Fiz. Nauk 
111, 3 (1975) [Sov. Phys. Usp. 18, 649 (1976)] 

[10] B. G. Englert, " Semiclassical Theory of Atoms". Lecture Notes in Physics, 
300 (Springer, 1988) 

[11] G. V. Shpatakovskaya, Teplofiz. Vysok. Temp. 23, 42 (1985) [Sov. Phys. 
High. Temp. 23, 36 (1985)] 

[12] E. A. Kuzmenkov and G.V. Shpatakovskaya, Int. J. Thermophys., 13, 315 
(1992) 

[13] V. Kresin, Phys.Rev. B38, 3741 (1988) 

[14] M. Membrado and A. F. Pacheco, Phys.Rev. B41, 5643 (1990) 

[15] V. M. Strutinsky, Nucl.Phys. A122, 1 (1968) 

[16] C. Yannouleas and Uzi Landman, Phys.Rev. B48, 8376 (1993) 

[17] D. A. Kirzhnits and G. V. Shpatakovskaya, Zh.Eksp.Teor.Fiz. 62, 2082 
(1972) [Sov.Phys.- JETP, 35, 1088 (1972)] 



6 



